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Abstract. Dark energy is the invisible fuel that seems to drive the current acceleration of the Universe. Its presence, which is 
inferred from an impressive convergence of high-quality observational results along with some sucessful theoretical predictions, 
is also supported by the current estimates of the age of the Universe from dating of local and high-z objects. In this work we 
study observational constraints on the dark energy equation of state (w) from lookback time measurements of high-z galaxies, 
as recently published by the Gemini Deep Deep Survey (GDDS). In order to build up our lookback time sample from these 
observations we use 8 high-z galaxies in the redshift interval 1.3 < z < 2.2 and assume the total expanding age of the Universe 
to be tg bs = 13.6^ 3 Gyr, as obtained from current cosmic microwave background data. We show that these age measurements 
are compatible with values of w close to -1, although there is still space for quintessence (w > -1) and phantom (w < -1) 
behaviors. In order to break possible degeneracies in the O m - w plane, we also discuss the bounds on this parametric space 
when GDDS lookback time measurements are combined with the most recent SNe la, CMB and LSS data. 
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1. Introduction 

Over the past decade, an increasing number of observational 
results have consistently indicated that we live in a nearly 
flat, accelerating universe composed of 1/4 of pressure- 
less matter (barionic + dark) and 3/4 of an exotic com- 
ponent with large negative pressure, usually called dark en- 
ergy or quintessence. The basic set of observations supporting 
such idea includes distance measurements to intermediary and 
high-z Type la Supernovae [SNe la] (Riess et al. 1998; 2004; 
Perlmutter et al. 1999; Astier et al. 2006), measurements of 
the Cosmic Microwave Background (CMB) anisotropies (de 
Bernardis et al. 2000; Spergel et al. 2003; 2006) and the cur- 
rent observations of the Large-Scale Structure (LSS) in the 
Universe (Tegmark et al. 2004; Eisenstein et al. 2005; Cole 
et al. 2005). Although it is widely believed that the existence 
of this dark energy component provides the key and remaining 
pieces of information connecting the inflationary flatness pre- 
diction (Q-Totai = 1) with astronomical observations, nothing but 
the fact that it has a negative equation of state (EoS) w = p/p 
and that its energy density and pressure are of the order of the 
critical density (i.e., p ^ -p ^ -10~ 29 g/cm 3 ) is known thus far 
(for recent reviews on this topic, see Peebles and Ratra, 2003; 
Padmanabhan 2004; Copeland et al. 2006). 
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On the other hand, it is also well true that there is now a 
flurry of activity focused on unveiling the nature of dark en- 
ergy or constraining any alternative explanation for the cosmic 
acceleration [e.g., extra dimensions effects (Sahni and Shtanov, 
2003; Zhu and Alcaniz, 2005)]. However, the absence of a nat- 
ural guidance from particle physics theories on the origin of 
this dark component seems to make clear that in order to bet- 
ter understand its actual nature, an important strategy is to find 
new observational methods or to revive old ones that could di- 
rectly or indirectly quantify the amount of dark energy present 
in the Universe, as well as determining its EoS. In this regard, 
the possibility of constraining cosmological parameters from 
age estimates of high-z objects constitutes an important and 
interesting attempt. In reality, since the days pre-dark energy, 
estimates of the age of the Universe (fu) have been one of the 
most pressing piece of data supporting such an idea (Krauss 
and Turner, 1995), in that dark energy helps explain the cur- 
rent dating of globular clusters by allowing a period of cosmic 
acceleration and leading to a larger expansion ageQ. In this re- 
gard, it also worth emphasizing that the evolution of the age 

1 For the widely accepted current value of the Hubble parameter, 
i.e., H a 72 km.s-'.lVqxr 1 (Freedman et al. 1998), no fiat CDM 
model without dark energy (whose age prediction is H B t\j = 2/3) may 
be compatible either with the direct age estimates from globular clus- 
ters or with the indirect age estimates, as provided by SNe la and CMB 
measurements (Tonry 2002; Spergel et al. 2003; 2006). 
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of the Universe with redshift (dtu/dz) differs from scenario to 
scenario, which means that models that are able to explain the 
total expanding age (at z = 0) may not be compatible with age 
estimates of high-z objects (and vice-versa). This in turn rein- 
forces the idea that dating of objects at high-redshift constitutes 
one of the most powerful methods for constraining the age of 
the Universe at different stages of its evolution (Dunlop et al. 
1996; Kennicutt Jr. 1996; Spinrad et al. 1997; Hasinger et al. 
2002), the first epoch of galaxy/quasar formation (Alcaniz & 
Lima 2001; Alcaniz et al. 2003), as well as for discriminating 
among different dark energy models (Krauss 1997; Alcaniz & 
Lima 1999; Jimenez & Loeb 2002; Jimenez et al. 2003; Fria^a 
et al. 2005; Jain and Dev, 2006). 

Recently, with the advent of large telescopes and new tech- 
nics, it was possible to estimate more precisely ages of high- 
z objects, including galaxies, quasars and galaxy clusters. In 
this regard, of particular interest to test current dark energy 
parametrizations are the discoveries of the prototypical evolved 
red galaxies LBDS 59W091 and LBDS 59W069 (Dunlop et al. 
1996; Spinrad et al. 1997; Bruzual and Magris 1997; Dunlop 
1999; Yi et al. 2000; Nolan et al. 2001), the quasar APM 
08279+5255 (Hasinger et al. 2002; Friaca et al. 2005) and the 
recently released sample of old passive galaxies from Gemini 
Deep Deep Survey [GDDS] (Abraham et al., 2004; McCarthy 
et al., 2004). In particular, the GDDS sample consists of 20 ob- 
jects lying in the interval 1.3 < z < 2.2, whose integrated light 
is dominated by evolved stars. As discussed by McCarthy et al. 
(2004), the data seems to indicate that the most likely star for- 
mation history is that of a single burst of duration less than 0. 1 
Gyr, although in some cases the duration of the burst is con- 
sistent with Gyr, which means that the galaxies have been 
evolving passively since their initial burst of star formation. 
From the theoretical viewpoint, these data, along with other age 
estimates of high-z objects, were recently used to reconstruct 
the shape and redshift evolution of the dark energy potential 
(Simon et al., 2005), as well as to place bounds on holography- 
inspired dark energy scenarios (Yi and Zhang, 2006). 

In the present paper, by following the methodology pre- 
sented by Capozziello et al. (2004) [see also Saini et al. 2000 
and Pires et al. 2006], we discuss quantitatively how GDDS 
age estimates of high-z galaxies constrain the EoS describing 
the dark energy. In order to perform our analysis, we transform 
the selected GDDS observations into lookback time (LT) mea- 
surements by assuming the total expanding age of the Universe 
to be t°]' s = 13.6^3 Gyr, as recently obtained from an analysis 
involving the most recent CMB data (MacTavish et al., 2005). 
To better constrain the dark energy EoS and energy density, we 
also investigate the bounds on these quantities when GDDS LT 
measurements are combined with the most recent SNe la and 
LSS data. 

2. Background Equations and Lookback Time 

With the assumption that the effective equation of state, 
w ~ f w(z)Qx(z)dz/£lx(z)dz, is a good approximation for a 
wide class of dark energy scenarios, we start our analysis by 
considering a homogeneous, isotropic, spatially flat universe 
described by the Friedmann-Robertson-Walker line element, 



ds 1 = dt 2 - a 2 (t)(dx 2 + dy 2 + dz 2 ), where a(t) is the cosmo- 
logical scalar factor and we have set the speed of light c = 1 . 
The Friedmann equation for such a scenario is given by 



•H(p) = VE(p) with E(p) = Q m a~ 3 + Q x a~ 3(1+w) 



(1) 



where "W(p) = H(p)/Ho, the complete set of parameters is 
p = {Q.j, w] (j = m and x stand for matter and dark energy 
density parameters, respectively), and from now on the sub- 
script denotes present-day quantities. In this background, the 
lookback time-redshift relation, defined as the difference be- 
tween the present age of the Universe (f ) and its age (f z ) when 
a particular light ray at redshift z was emitted, can be written as 



;p) = h - 1 r 

Jo 



dz' 



(l+zOW 



(2) 



where H ' = 9.78/z 1 Gyr and h ranges in the interval 0.64 < 
h < 0.8, as provided by the HST key project (Friedman et al., 
2002). 

To proceed further, let us now consider an object (e.g., a 
galaxy, a quasar or a galaxy cluster) at redshift z, whose the 
age t(zi) is defined as the difference between the age of the 
Universe at z, and the one when the object was born (at its 
formation redshift zf ), i.e., 



tizd - k 



C°° dz' f°° dz!_ 

J Zi (1+z'Wp) J ZF d+z'Y. 

or, equivalently, 



)W(p) 



t(zd = t L {z F ) - Hzd. 



(3) 



(4) 



From the above expressions, we can define the observed look- 
back time to an object at z, as 



ff 5 fc;r) = t L {z F )-t{ Zi ) 
■ Kzdl 



r f obs 
L'o 

+obs 



L'o 



Hzf)] 



■ KZi) - T, 



(5) 



where t stands for the incubation time or delay factor, which 
accounts for our ignorance about the amount of time since the 
beginning of the structure formation in the Universe until the 
formation time (//) of the object [see also Capozziello et al. 
(2005) and Pires et al. (2006) for a similar discussion]. 

3. Lookback time Data and Statisical analysis 

3.1. Data 

In this section we discuss quantitatively how the GDDS LT 
estimates of high-z galaxies constrain the EoS describing the 
dark energy through a statistical analysis of the data. The to- 
tal GDDS sample presented by McCarthy et al. (2004) consists 
of 20 old passive galaxies distributed over the redshift inter- 
val 1.308 < z < 2.147 and reveal unambiguous evidence for 
old and metal-rich galaxies over the full redshift range. To de- 
termine the galaxy ages McCarthy et al. (2004) compare the 
observed spectral energy distribution with a set of synthetic 
spectra computed with PEGASE.2 (Fioc & Rocca-Volmerange, 
1997) and construct a multidimensional^ 2 surface spanning a 
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Fig. 1. Left: The age-redshift data points, (a) Original data from GDDS. This sample corresponds to 20 old passive galaxies 
distributed over the redshift interval 1.308 < z < 2.147, as given by McCarthy et al. (2004). (b) The 8 high-z measurements 
selected after the criterion discussed in the text. Right: Lookback time relation as a function of the redshift for some values of 
the EoS parameter w and a fixed value of Q m = 0.27. To plot these curves, we have used r = 1.5 Gyr. Thick line corresponds to 
the ACDM case (w = -1). 



wide range of star formation histories, extinction and metal- 
icities. Although this galaxy sample has been independently 
reanalyzed by Simon et al. (2005), who obtained ages within 
0. 1 Gyr of the GDDS collaboration estimates, it is important 
to emphasize that these age estimates are not free of observa- 
tional and theoretical uncertainties (for instance, those related 
to the population synthesis models, the amount of active evolu- 
tion, among others.) [As an example, see the debate involving 
the age estimates for the radio galaxy LBDS 59W069; Bruzual 
and Magris, 1997; Yi et al, 2000; Nolan, Dunlop and Jimenez, 
2001]. 



In order to build up our LT sample, we first select from 
GDDS observations what we will consider as the most appro- 
priate data to our cosmological analysis. Since our primary ob- 
jective here is to impose limits on the behavior of the cosmic 
EoS as restrictive as possible from these high-z age measure- 
ments, we adopt the criterion that given two objects at (approx- 
imately) the same z, the oldest one is always selected. This is 
justified by the fact that a universe model that is capable of 
explaining the existence of a very old object at a given z can 
also naturally explain the age of the youngest objects at that z. 
By following this criterion we end up with a sample of 8 data 
points, which are shown in Fig. 1(b) along with the original 
data from GDDS collaboration [Fig. 1(a)]. The other important 
aspect to build up our lookback time sample [see Eq. (5) above] 
concerns the total age of the Universe, which is assumed in our 
analysis to be f° 6s = 13. Gyr, as obtained by MacTavish 
et al. (2005) from a joint analysis involving only data of the 



most recent CMB experiments (WMAP, DASI, VSA, ACBAR, 
MAXIMA, CBI and BOOMERANG) H 

3.2. Statistical Analysis 

In order to estimate the best-fit to the set of parameters p we 
define the likelihood function 

tage «: exp [-^fo p, r)/2] , (6) 
where xl ge is given by 

2 _ v [fcfc;p)- f L bs (zr, T)f 

. u T 

[Up) - t° hS f 
+ - 2 L - (7) 

°%> 

Here, (Tj = crj + <r^, bs , <rj is the uncertainty in the individual 

lookback time to the i th galaxy of our sample and cr t „i,.< - 0.35 
Gyr stands for the uncertainty on the total expanding age of the 
Universe (f bs ). As mentioned earlier, the evolution of the age 
of the Universe with redshift may differ considerably from sce- 
nario to scenario, which amounts to saying that it is possible 
that cosmological models that are able to explain age estimates 
of high-z objects may not be compatible with the total expand- 
ing age at z = (Friaca et al. 2005). In our statistical analysis, 

2 Note that, to avoid double counting of information with SNe la 
and LSS analyses performed in the next sections, we have adopted 
age estimates obtained from combinations of CMB experiments only. 
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Fig. 2. First results of our statistical analyses. GDDS LT con- 
straints on the Q m - w plane. The contours correspond to 68.3% 
and 95.4% confidence levels. The best-fit model for this analy- 
sis are £2 m = 0.21 and w = -0.72. 

this is taken into account by including a prior on the total age of 
the Universe, which is equal (as should be) to the value adopted 
to build up the lookback time sample, i.e, f° 6j = 13.6*q1 Gyr. 
This in turn amounts to saying that our LT analysis is always 
performed with the two terms of Eq. ||7). 

Another important aspect concerns the delay factor t. Note 
that in principle there must be variations in the value of r for 
each object in the sample (galaxies form at different epochs). 
Here, however, since we do not know the formation redshift 
for each object, the delay factor is assumed as a "nuisance" 
parameter, so that we marginalize over it. Note also that this 
marginalization may also be analytically obtained by defining 
a modified log-likelihood function^ 2 , i.e., 



x 1 



-2 In 



f 

Jo 



drexpl-i^; 



(8) 



B 2 

A + £> - 2 In 

C 



£ erfc (\jy 



where 



rr 2 ' ~ ^ a 2 ' ~ ^ cr 2 ' 

D is the second term of the rhs of Eq. ©, 
A = fefc;p) -[/?'-/&)], 

and erfc(jc) is the complementary error function of the variable 
x. 

4. Results 

Figure l(c) shows the GDDS LT data as a function of the red- 
shift for some selected values of the EoS parameter and a fixed 



value of Q m = 0.27. To plot these curves we have assumed a 
delay factor of t = l .5 Gyr. Note that, in agreement with other 
independent analyses, only for extreme (and unrealistic) val- 
ues of the delay factor (e.g., t > 4.5 Gyr), a matter-dominated 
universe (w = 0) would be still compatible with the GDDS LT 
estimates. 

In Fig. 2 we show the first results of our statistical analysis. 
Confidence regions (68.3% and 95.4% c.l.) are shown in the 
parametric space Q. m - w from the GDDS LT data discussed 
above. Similarly to the present constraints from SNe la (see, 
e.g., Astier et al., 2006), we note that a large interval for the 
dark energy EoS is allowed, whereas at 2cr level the matter den- 
sity parameter is more restricted here (Q m < 0.35) than in the 
current SNe la analyses (Q. m < 0.5). Note also that, although 
compatible with a supernegative behavior of the dark energy 
EoS (w < -l), the best-fit value for w is much less negative 
(w = -0.72) than the one predicted by the so-called phantom 
scenarios (see, e.g., Caldwell, 2002; Faraoni, 2002; Alcaniz, 
2004 Feng, Wang and Zhang, 2005 and references therein). For 
the sake of comparison, we add that If we had used the total 
GDDS sample of 20 galaxies (instead of applying and using 
the criteria discussed earlier to reduce the number of objects), 
our best-fit result would be Q. m =s and a> - -0.32, which 
is very far from any realistic estimate using the current sets of 
observational data. It is also worth observing that GDDS LT 
data provide contours in the Q. m - w plane very similar to those 
arising from distance meausurements from SNe la observations 
(see Fig. 3b), which means that a joint analysis involving these 
data sets may also be understood as an extended one, in the 
sense that while SNe la data reach z < 1 .01 (Astier et al., 2006), 
GDDS data lies in the redshift interval 1.308 < z < 2.147. 

4. 1 . Complementary Analysis 

In order to complement our lookback time study on the para- 
metric plane O m - w, we combine our GDDS LT data with 
some of the most recent sets of cosmological observations. 
We use here the current SNe la data from Supernova Legacy 
Survey (SNLS) collaboration corresponding to the first year re- 
sults of its planned five years survey (Astier et al. 2006). The 
SNLS sample includes 71 high-z SNe la in the redshift range 
0.2 < z < 1 and 44 low-z SNe la compiled from the literature 
but analyzed in the same manner as the high-z sample. This 
data-set is arguably the best high-z SNe la compilation to date 
and seems to be in good agreement with WMAP results in that 
it favors values of w ^ -1 (Jassal et al. 2006). We also include 
in our joint analysis current LSS data to help break the degener- 
acy between the dark energy EoS and the matter density param- 
eter. For the LSS data, we use the recent measurements of the 
BAO peak in the large scale correlation function detected by 
Eisenstein et al. (2006) using a large sample of luminous red 
galaxies from the SDSS Main Sample. The SDSS BAO mea- 
surement provides m = 0.469(« s /0.98)" 35 + 0.017, with 31 
defined as 



1/2 



-ZBAO 



Zbao 



r 2 (z B Ao;p) 



E(z B ao;p) 



1/3 



(9) 
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Fig. 3. Complementary results of our statistical analyses, a) GDDS LT plus SNLS constraints on the Q m - w plane. The contours 
correspond to 68.3% and 95.4% confidence levels. The best-fit model for this analysis are Q m = 0.28 and w = -1.05 with 
X 2 /v = 1.02 . The dotted lines represent the constraints from BAO measurements, b) SNLS constraints on the parametric space 
f2 m - w shown here for the sake of comparison (see also Astier et al. 2006). c) The same as in Panel (a) when GDDS LT data are 
combined with BAO measurements, d) Current constraints on the Q m - w plane from SNLS + BAO data. 



where T is the distance coordinate, zbao = 0.35, E(zbao; p) is 
given by Eq. (1), and we take the scalar spectral index n$ = 
0.95, as given by Spergel et al. (2006). 

Figure 3(a) shows the resulting likelihood from the GDDS 
LT + SNLS data. There, confidence contours at 68.3% and 
95.4% are shown in the plane Q. m - w. Note that the allowed 
parametric space is now considerably reduced relative to the 
previous panel (Fig. 2) and to Fig. 3(b) which shows, for 
the sake of comparison, the Q. m - w constraints from SNe la 
alone (see also Astier et al. 2006). In particular, less negative 
(w > -0.7) and very negative (w < -1.75) values of w are ex- 
cluded at 2<x level. Note also that while the SNLS data alone 
cannot place any restrictive bounds on the dark energy EoS at 
1 or 2cr levels, the combination GDDS LT + SNLS data moves 
up the confidence regions to values of w - -1, which is in 
agreement with other current sets of data (see, e.g, Jassal et al. 
2006 for a discussion). 

The dotted lines in Fig. 3(a) represent, respectively, the con- 
straints from BAO measurements on the plane £2 m - w. As is 
well known, since this quantity has been measured at a spe- 
cific redshift (zbao = 0.35), it forms bands on the parametric 
space Q. m - w instead of ellipsoids, as in the case of GDSS 
lookback time and SNe la data. Here, however, the most im- 
portant aspect associated with these measurements lies in the 



fact that the bands of allowed parameters are roughly perpen- 
dicular to the major axis of the ellipsoid arising from the GDDS 
LT + SNLS data, which suggests that a combined analysis may 
go even further in breaking the degeneracies inherent to the 
parametric plane O m - w. The results of these complementary 
analyses are shown in Figs. 3(c), 3(d) and 4. In the former, we 
show the GDDS LT + BAO constraints on the plane Q m - w. 
Note that a considerable part of the parametric space is now 
ruled out at 2<x level [mainly relative to Panel 3(a)] and an up- 
per (lower) limit w < -0.75 (w > -1.4) can be roughly placed 
on the EoS parameter. These results should be compared with 
the very restrictive constraints from SNLS + BAO data [Panel 
3(d)], providing Q m = 0.271 and w - -1.023 (Astier et al. 
2006). As expected, when all the cosmological observables dis- 
cussed in this paper are combined in a joint analysis (Fig. 4) 
the allowed parametric space is even more reduced relative to 
the previous analyses (Figs. 2, 3a-3d), with the value of dark 
energy EoS converging to ^ - 1 (although there is still space 
for quintessence and phantom behaviours). The best-fit param- 
eters for this total joint analysis are O m = 0.27 and w — -1.02 
with reduced^ = 1.02 (= X 2 m iJ v ^ where v is the numeber 
of degrees of freedom). This roughly corresponds to the cur- 
rent standard ACDM scenario, i.e., an accelerating universe 
with deceleration parameter qo = -0.57, transition redshift (at 
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Fig. 4. Joint analysis involving GDDS LT + SNe la + BAO 
data. Note that, relative to the previous Panels, the allowed 
parametric space is now considerably reduced. The best-fit val- 
ues for all the analyses performed in this paper are presented in 
Table!. 



Table 1. Best-fit values for Cl m and w 



Test 


Q. m 


w 


X 2 lv 


GDDS Lookback Time 


0.21 


-0.72 


1.55 


GDDS Lookback Time + SNLS 


0.28 


-1.05 


1.02 


GDDS Lookback Time + BAO 


0.29 


-1.03 


1.58 


Lookback Time + SNLS + BAO 


0.27 


-1.02 


1.02 



which the expansion switches from deceleration to accelera- 
tion) zj = 0.75, and total expanding age of 9.6/z -1 Gyr 13.5 
Gyrfor/z = 0.71). At 95.4% c.l. we also found Q m = 0.27+0.02 
and w = -1.02+qJj. The best-fit parameters for all the analyses 
performed in this paper along with the corresponding xl are 
summarized in Table I. 

5. Summary 

The recent accumulation of independent observational results 
has opened up a robust window for probing the behavior of the 
dark energy. However, as is well known, most of the methods 
employed to place limits on the dark energy EoS (w) or, more 
generically, on the parametric space Q m - w, are essentially 
based on distance measurements to a particular class of objects 
or physical rulers (SNe la, CMB, galaxy clusters, etc.). In this 
regard, it is also particularly important to obtain accurate and 
independent bounds on the physical behavior of the dark en- 
ergy, as well as on the other main cosmological parameters, 
from physics relying on different kinds of observations. 

In this paper we have followed this direction and studied 
constraints on the parametric space O m - w from age measure- 
ments of high-z galaxies, as recently provided by the Gemini 



Deep Deep Survey. By transforming the GDSS data into look- 
back time estimates (by using current estimates of the total age 
of the Universe from LSS and CMB data), we have discussed 
quantitatively how the current age data constrain the EoS de- 
scribing the dark energy. We have shown that, although allow- 
ing a large interval for w (similarly to current SNe la data), 
these LT estimates prefer values of w > -1 (best-fit) and con- 
strain the matter density parameter to be £2 m < 0.35 at 2cr level. 
The most restrictive bounds on the parametric plane £2 m - w 
are obtained when age and distance measurements are com- 
bined in a joint statitical analysis. In this case, using GDDS LT 
along with the most recent SNe la and LSS observations, we 
have found Q m = 0.27 ± 0.02 and w = at 95.4% 

(c.l.) which, although being compatible with a quintessence 
and phantom cases, clearly favors the cosmological constant 
behavior. 

Finally, we emphasize the importance of more and more 
precise age measurements of high-z objects. These observa- 
tions will certainly provide a new and complemetary tool to test 
the reality of the current cosmic acceleration, to place restric- 
tive bounds on the cosmic EoS characterizing the dark energy 
and, as a consequence, to distingush among the many alterna- 
tive world models. The present work, therefore, highlights the 
cosmological interest in the observational search for old col- 
lapsed objects at low, intermediary and high redshifts. 
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